Three Dimensional Magnetism and Coupling to the Conduction Electrons in PdCr02 
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We report density functional calculations addressing the electronic structure and magnetic prop- 
erties of delafossite PdCr02. We find substantial magnetic interactions in the c-axis direction as 
well as beyond first nearest neighbors in-plane, so that PdCr02 is a 3D frustrated antiferromagnet. 
We also find substantial coupling between the Cr moments and the Pd derived conduction electrons. 

PACS numbers: 75.10.Lp,75.30.Et 



I. INTRODUCTION 

Recent discoveries pointing to novel ground states and 
unusual behaviors in field and temperature dependent 
magnetic and thermodynamic properties have led to re- 
newed interest in the physics of frustrated spin sys- 
tems. Delafossite structure PdCr02 has attracted re- 
cent attention in this context. Importantly, like the sis- 
ter compound PdCo02, it can be grown in high quality 
form facilitating experimental investigation of its phys- 
ical properties. This material can be described as a 
stacking of 2D Cr02 layers, separated by Pd, with the Cr 
ions arranged on a triangular lattice. In a standard near- 
est neighbor Heisenberg picture, the ground state would 
be expected to be a non-collinear 120° spin structure be- 
low the Neel temperature, TAr=37.5 K,^ but the actual 
magnetic structure of PdCr02 could be more complex.^ 
Takatsu and co-workers find that the magnetic Bragg 
peaks in neutron scattering are rather broad and half- 
integer / (1/3,1/3,/) and (2/3,2/3,/) peaks are found in 
addition to the integer / peaks. ^ In any case, PdCr02 
is a particularly interesting system because like its sis- 
ter compound PdCo02 it is conducting^ and has a very 
large conductivity anisotropy.^ This offers an opportu- 
nity to study the interplay between frustrated magnetism 
and charge carriers in a near two dimensional metal. Fur- 
thermore, recent experiments have revealed an unconven- 
tional anomalous Hall effect in this material.^ 

Doumerc and co-workers reported susceptibility data 
up to ~ 450 K.^ They obtained an effective moment of 
4.1 /iB and a Weiss temperature of -500 K, in a Curie- 
Weiss fit, although they qualify these values as approxi- 
mate. The data show a poor fit except over a relatively 
small temperature range above ~ 300 K, in contrast to 
the other compounds they report: CuCr02, AgCr02 and 
CuFe02, which show very good Curie- Weiss behavior 
over a wide temperature range extending to below 100 
K. In fact, in their data, there is still an apparent up- 
ward curvature to x~^(^) at the highest T, suggesting 
that the Weiss temperature could be lower in magnitude 
or that the system may have itinerant character. 

Takatsu and co-workers show more details of the sus- 



ceptibility below 350 K.^ There is no plausible Curie- 
Weiss fit to the data, but they do state that their data 
are consistent with the values extracted by Doumerc and 
co-workers. Interestingly, the data shows a broad max- 
imum above T/v- However, in contrast to strongly 2D 
magnetic systems, x ji^st above Tn is only slightly lower 
than the maximum value, and there are strong signatures 
of the ordering in x(T) as well as resistivity and specific 
heat. This is not really consistent with expectations for 
a strongly 2D Heisenberg system. On the other hand, 
Takatsu and co-workers do observe a sublinear resistiv- 
ity, p{T) above Tat, which they attribute to short range 
magnetic correlations, and an unusual behavior in the 
specific heat also above Tn- 

Motivated by this we performed density functional 
calculations of electronic structure and energetics of 
PdCr02 in order to address the nature and role of the 
metallic conduction electrons in PdCr02 and the extent 
to which the system is a realization of a 2D nearest neigh- 
bor Heisenberg magnet. We find in contrast to previous 
assumptions that (1) there is a considerable interplay be- 
tween the Cr moments and the metallic electrons, and 
(2) there are strong magnetic interactions along the c- 
axis direction, so that even though PdCr02 is very two 
dimensional as an electron gas, it is a very three dimen- 
sional magnetic system. Finally, there is metal-metal 
bonding, consistent with previous studies of delafossite 
compounds 



II. APPROACH 

We performed density functional calculations with the 
Perdew, Burke, Ernzerhof (PBE) generalized gradient 
approximation,^^ using the all electron linearized aug- 
mented planewave (LAPW) method^ ^ as implemented 
in the WIEN2k code,^^ similar to our previous calcula- 
tions for PdCo02.^^ We carefully tested the convergence 
of our results against the various parameters, including 
tests with different sphere radii and different choices of 
the augmentation including both LAPW and so-called 
APW+lo basis sets.^^ We found the results to be stable. 
We used the experimental lattice parameters, a=2.923 A 
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FIG. 1. (color online) Density of states and projections on Cr 
and Pd LAPW spheres of radius 2.05 bohr for ferromagnetic 
ordering 

and c=18.087 A (delafossite structure, space group R3m, 
Pd at (|,^,^), Cr at (0,0,0), O at {±z,±z,±z), in rhom- 
bohedral coordinates).^^ Refinement of the internal pa- 
rameter corresponding to the O height above the Cr plane 
has not been reported in literature to our knowledge. As 
such, we calculated this parameter for a ferromagnetic 
ordering and used it in the other calculations. We ob- 
tained z=0.6111, which yields a Cr-0 distance of 1.962 
A. This is in accord with the expected bond length for 
high spin Cr^+ (the Shannon radii are 1.22 A for O^- and 
0.755 A for octahedral Cr^+, summing to 1.975 A). We 
also note that the Pd-Pd distance of a=2.923 A is not 
much longer than the Pd-Pd distance of 2.75 A in Pd 
metal. The Cr-O-Cr bond angles in the layers are 96.3°. 
Considering the Cr lattice, the O atoms directly connect 
the nearest neighbor Cr atoms in plane. There are no di- 
rect O connections to second neighbors. Secondly, in the 
out-of-plane direction, the shortest hopping paths are Cr- 
0-Pd-O-Cr. The Cr-O-Pd bond angles along this path 
are 120.7° with this O position. 



III. MOMENT FORMATION AND 
ELECTRONIC STRUCTURE 

We find that PdCr02 is strongly unstable against Cr 
moment formation regardless of the arrangement of the 
moments with an energy of ~1 eV/Cr. We also find 
all arrangements are metallic in accord with experiment. 
Therefore, PdCr02 should be described as a metal con- 
taining Cr local moments. We start our description with 
the calculated electronic structure with ferromagnetic 
alignment of these moments. The calculated energy is 
0.92 eV lower than the non-spin-polarized (no Cr mo- 
ments) case, and the magnetization is 2.87 /j^b per for- 
mula unit. The magnetization inside a Cr LAPW sphere 
(radius 2.05 bohr) is 2.57 /j^b- 

The calculated electronic density of states (DOS) is 
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FIG. 2. (color online) Fermi surface for ferromagnetic order- 
ing. Left panel is majority spin, right panel is minority. Note 
that these are hole cylinders centered at F. 



given in Fig. [H It clearly shows crystal field split Cr d 
states with the majority spin t2g states occupied, and the 
Cg unoccupied, as are all the minority spin Cr states. The 
DOS at Ep is Pd derived. These states provide a simple 
very two dimensional Fermi surface (Fig. [2j), consistent 
with the observed strong conductivity anisotropy and 
reminiscent of PdCo02, which is also a highly anisotropic 
metal. As discussed previously, the states giving rise to 
the Fermi surface arise from Pd-Pd metal-metal bonding 
as in other delafossites.*^ Additionally, two additional 
features are seen: (1) The Fermi surface has very little 
rounding, which is the condition for maximizing nesting 
(but note that density of states is Pd derived and low, 
N{Ef)=0.7 eV~^, per formula unit on a both spins basis 
for this ferromagnetic ordering), and (2) While both the 
majority and minority spin Fermi surfaces are similar in 
shape, they are not identical, implying an interplay be- 
tween the metallic, largely Pd derived, conduction band, 
and the Cr derived magnetic moments. In fact, the mi- 
nority spin Fermi surface is 13% larger in volume than 
the majority spin surface. This interplay between Cr 
moments and the conduction electrons provides a mech- 
anism for spin scattering, that will increase the resistivity 
as T disorders the spins and also a mechanism for long 
range Cr-Cr interactions in-plane, as well as interactions 
in the c-axis direction. 



IV. MAGNETIC INTERACTIONS 

We studied the magnetic interactions through a series 
of supercell calculations with different magnetic order- 
ings. These calculations and their energetics are sum- 
marized in Fig. [3l One may immediately note from the 
top panel of the figure that the out-of-plane magnetic 
interactions are large. Specifically, an antiferromagnetic 
stacking of ferromagnetic layers is found to be 18 meV/ Cr 
lower in energy than a ferromagnetic stacking. This is 
large compared with the ordering temperature, Tn= 38 
K, i.e. kTN= 3.3 meV. 

Of the small collinear cells that we studied, the low- 
est energy belongs to AF-2 (see Fig. [3] for notation), 
which corresponds to the ordering in CuFe02J^ — The 
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FIG. 3. (color online) Structures and energetics of small su- 
percells, in meV per Cr, relative to the ferromagnetic or- 
dered configuration. The top panel is for ferromagnetic layers, 
stacked along the c-axis either ferromagnetically, or antifer- 
romagnetically, while the bottom panel shows three in-plane 
antiferromagnetic orderings. The unit cells are indicated by 
dotted black lines. The three numbers in parentheses below 
each structure denote the numbers of opposite spin neighbors, 
at the the nearest neighbor position (out of six), the in-plane 
next nearest Cr neighbor position (out of six) and the out-of- 
plane nearest neighbor position (out of six), respectively. 
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FIG. 4. (color online) Fermi surface viewed down the c-axis 
for the lowest energy AF-2 configuration (see text), shown 
in an extended zone scheme. The white lines show the su- 
percell Brillouin zone and the highlight sketches the non- 
reconstructed Fermi surfaces of Fig. [2l 



electronic structure for the various the orderings remains 
very anisotropic at the Fermi energy, where the states 
are Pd derived. However, as mentioned above, there is 
a strong interplay between the Fermi surface and the Cr 
moments. For example, the Fermi surface for the lowest 
energy magnetic configuration is shown in Fig. |4j The 
figure shows overlapped hexagonal sections, as expected 
for the folded zone, but also sizable splittings at the in- 
tersections. 



This interplay between magnetic order and the con- 
duction electrons is also seen in the calculated densities 
of states. As mentioned, for the ferromagnetic order, we 
obtain N{Ef)=0.70 eV~^ on a per formula unit both 
spins basis. For the AF-z (see Fig. [3j) order, we obtain 
N{Ef)=0.69 eV~^ on the same basis, while we obtain 
N{Ef)= 0.81 eV-\ for AF-1, N{Ef)= 0.76 eV-\ for 
AF-2, and N{Ef)= 0.77 eV-\ for AF-3. The bare spe- 
cific heat coefficient 7 inferred from these values is 1.6- 
1.9 mJ/mol K^, comparable to the specific heat value 
of Takatsu and co-workers of 1.4±0.2 mJ/mol in the 
ground state. ^ This leaves little room for any renormal- 
ization at low T. 

The variability of the electronic structure dit Ef de- 
pending on magnetic state provides a qualitative frame- 
work for understanding why the resistivity of high qual- 
ity crystals of PdCr02 is so much higher than that of 
PdCo02, even though the Fermi surface, structure and 
band character at the Fermi surface are very similar, ^^li^ 
and also why the resistivity shows a strong signature of 
the magnetic ordering, with a pronounced decrease in re- 
sistance as T is lowered through T^} Specifically, this 
connection between magnetic order and the electrons at 
the Fermi surface indicates that strong spin-fluctuation 
scattering is expected above the ordering temperature, 
freezing out as T is lowered below T^. 

Materials with this type of coupling between mag- 
netism and electrons at the Fermi surface tend to be the 
ones that display unusual properties, such as supercon- 
ductivity (as in e.g. the Fe-pnictides where the coupling 
is very strong) when the magnetism is suppressed. There- 
fore, this coupling of spin fluctuations to electrons at the 
Fermi surface suggests that it will be of considerable in- 
terest to experimentally examine what happens as the 
magnetic order is suppressed, e.g. by pressure or alloy- 
ing. 

We now turn to the energetics (Fig. [3]) in more detail. 
To do this we write a short range model with interactions 
to nearest and next-nearest neighbors in-plane and to the 
nearest out-of-plane neighbor: 

E = E^F) + jiNn + j2N2i + jzN,^ (1) 

where E is the energy per Cr, ^(f) is the energy of the 
ferromagnetic ordering, is the average number of op- 
posite spin first nearest neighbors in the structure, and 
similarly for next nearest neighbors and neighbors along 
the c-axis direction. While this could be converted into 
a three neighbor Heisenberg model, we write it this way 
because the moments in the Cr LAPW spheres vary be- 
tween the different configurations (in the range 2.48 fiB 
to 2.57 ijlb)' 

In any case, if we match the energies of the ordered 
structures, F, AF-z, AF-1 and AF-3, we obtain ji= -16 
meV, j2= +2 meV and jz= -3 meV. The energies cannot 
be reasonably reproduced without all three parameters 
(note that each shell has 6 atoms, so e.g. the difference 
between ferromagnetic and AF-z is 6^2=- 18 meV, which 
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is not small, etc.) Use of these values, yields the energy 
for the AF-2 state to within 1 meV / Cr atom. While 
the model is perhaps too simple to capture all aspects 
of the magnetic interactions in PdCr02, it does clearly 
show that the material cannot be reasonably understood 
as a simple 2D nearest neighbor Heisenberg model. 



V. DISCUSSION AND CONCLUSIONS 

Considering the ordering temperature /cT/v=3.3 meV, 
it would seem that none of these interactions can be ne- 
glected, i.e. at least first and second in-plane neighbors 
and crucially c-axis interactions need to be included when 
discussing the magnetism of PdCr02. 

This conclusion may seem unexpected considering the 
very two dimensional electronic structure at the Fermi 
energy, as found both in our calculations and as is clear 
from the experimental resistivity anisotropy.^ However, 
we note that a similar behavior is found in the layered 
cobaltate Na^ECo02.^^'^^ That materials has a related 
crystal structure, that is similar to the delafossite, ex- 
cept that the layer stacking differs and bridges between 
planes are O neighbors along c instead of two O sepa- 
rated by a Pd (i.e. in Naa;Co02 the O in neighboring 



layers are directly on top of each other). The three di- 
mensional magnetic behavior in Naa;Co02 arises in part 
because of this bonding topology, which provides many 
superexchange exchange paths between the neighboring 
Cdl2 structure C0O2 layers. Within this framework, 
the c-axis magnetic interactions and the coupling to the 
conduction electrons are inter-related since they are both 
mediated by or through Pd. 

In any case, we note that the Kosterlitz-Thouless type 
of suppression of the ordering temperature is logarithmic 
in the anisotropy of the interactions. As such, very strong 
anisotropy is needed to obtain two dimensional, as op- 
posed to three dimensional, magnetism. We can conclude 
that PdCr02 is a three dimensional frustrated antiferro- 
magnet with an interesting interplay between magnetism 
and conduction electrons. 
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